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ABSTRACT 



^^ , This review presents a brief discussion of the theory of Cold Dark Matter of 

_. I structure formation in the universe, describing the main processes determining the 

CSJ I power spectrum, the non-linear gravitational collapse, the formation of galaxies and 

the evolution of the intergalactic medium. Recent advances in the observations of 

^ ■ high redshift galaxies, and their interpretation in the context of this theory, are then 

2_? • summarized. 
l/^ . Este articulo de revision presenta una breve exposicion de la teoria de la Materia 



Invisible Fria de la formacion de estructura en el universe, describiendo los procesos 



Q^ I mas importantes que determinan el espectro de potencia, el colapso gravitatorio no 

5^ I lineal, la formacion de galaxias y la evolucion del medio intergalactico. Resumimos 

i-^ I tambien avances recientes en las observaciones de galaxias a alto corrimiento al rojo, y 

I ' la interpretacion en el contexto de esa teoria. 

^ ■ Subject headings: 

> 

^- 

j^ ■ 1. Introduccion 

Ciertas cuestiones que han resultado en grandes avances en la historia de la astronomia y la 
cosmologia tienen, ademas, un interes especial para la humanidad porque tratan sobre fenomenos 
totalmente imprescindibles para nuestra existencia en el universo: la fuente de energia del Sol, la 
formacion de sistemas planetarios, y la sintesis de los elementos en el interior de las estrellas son 
ejemplos que acuden a la mente. La cosmologia moderna, despues de habernos encaminado por 
vez primera al estudio cientifico del origen del universo observable en su totalidad, en el marco 
de la teoria de la Gran Explosion, trae consigo tambien otra cuestion dentro de esa categoria: 
el origen de las fluctuaciones primordiales y la formacion de las galaxias. El universo debe 
empezar con condiciones iniciales practicamente homogeneas para resultar en el estado presente 
de homogeneidad a gran escala; pero la existencia de fluctuaciones de densidad primordiales, de 
una amplitud 6p/p -^ 10~^, es absolutamente necesaria para la formacion de las galaxias. La 
gravedad puede amplificar y llevar las fluctuaciones iniciales al colapso no lineal, pero no puede 



crearlas. El proceso que genero esas fluctuaciones en el universo primitivo, sumido todavia en el 
misterio, dio lugar a la gran diversidad del universo no lineal, incluida la existencia de la vida, y 
evito la continuacion indefinida del universo lineal y homogeneo, conteniendo linicamente un mar 
de radiacion y atomos de hidrogeno y helio. 

Estamos actualmente en la era de exploracion y descubrimiento hacia la frontera de alto 
corrimiento al rojo o, en otras palabras, de las mayores distancias des de las que es posible 
recibir mensajes en el universo. Dada la velocidad de la luz, las grandes distancias proporcionan 
a la cosmologia la evidencia de la historia pasada de la formacion de las galaxias. Durante las 
ultimas tres decadas, nuestra vision del pasado del universo ha progresado continuamente con el 
descubrimiento de galaxias y micleos activos a distancias cada vez mayores, el estudio de fondos 
de radiacion cosmica a distintas frecuencias, y el analisis de espectros de absorcion del hidrogeno 
intergalactico interpuesto en la direccion de fuentes luminosas. Presentamos en este articulo de 
revision un breve resumen del estado actual de la teoria y observaciones sobre la formacion de las 
galaxias y la evolucion del medio intergalactico. No podemos hacer justicia en ese breve articulo 
al inmenso mimero de trabajos publicados en diversos temas de gran impacto sobre el estudio 
de la formacion de las galaxias. Entre otros articulos de revision, Ellis (1997) describe el estudio 
de galaxias debiles, y Ranch (1998) presenta una excelente exposicion de observaciones y teorias 
del bosque Lya de hidrogeno intergalactico. Cabe recomendar asimismo varios articulos en libros 
de conferencias recientes de gran utilidad para ponerse al dia en ese campo: vease, por ejemplo, 
Madau 1999, Steidel 1998a,b. 



2. Teoria de la Materia Invisible Pria 

De entre los distintos modelos propuestos de formacion de galaxias, la teoria de la Materia 
Invisible Fria (que abreviaremos MIF; en ingles. Cold Dark Matter) ha resultado ser la de mayor 
exito, y claramente favorecida por las observaciones. La teoria postula que la materia invisible, 
cuya presencia se deduce de las determinaciones dinamicas de la masa de galaxias y ciimulos (p.e., 
Trimble 1987), consiste en objetos o particulas "frias", sin dispersion de velocidades inicial, y que 
las fluctuaciones primordiales de densidad son adiabaticas (es decir, manteniendose constante 
la razon de la densidad de fotones, bariones y materia invisible), Gausianas, y con un espectro 
de potencias invariante de escala (vease por ejemplo, Blumenthal et al. 1984, Ostriker 1993). 
Generalmente, se supone que las fluctuaciones fueron generadas en un periodo de inflacion por 
algiin proceso que se mantuvo constante mientras el rango de escalas observable en el presente 
cruzaba el horizonte de acontecimientos, lo cual implica la invariancia de escala. 

Inicialmente, la teoria MIF se considero de forma casi exclusiva dentro del modelo cosmologico 
con densidad de materia critica, Q = p/pcrit = 1 (donde p es la densidad media de materia y 
pcrit es la densidad critica). Desde los inicios de la cosmologia observacional, las medidas de 
la densidad media del universo indicaron un valor Q < 1 (vease, p.e. Gott et al. 1974). Este 
resultado observacional se ha mantenido hasta la actualidad: la densidad media de luz en la banda 



B es en el presente pB — 2 x 10^ HLq Mpc^ (Zucca et al. 1997 y referencias incluidas; utilizamos 
aqui Hq = lOO/ikms^ Mpc^ ), y la razon masa-luminosidad en diversos sistenias colapsados, 
desde galaxias individuales a los niayores ciimulos, niedida a radios suficientemente grandes, tiende 
generalmente a valores 150/i ^ M /Lb ^ 500/i (p.e., Bahcall, Lubin, & Dorman 1995), resultando 
en una densidad de materia p ~ 5 x lO^'^/i^M© Mpc~^, y pj pent = i^ — 0.2. 

El valor de fi asi deducido a partir de la densidad luminosa del universo puede verse afectado 
por el efecto del sesgo de las galaxias. En principio, la mayor parte de la niasa del universo podria 
ubicarse en los grandes vacios, sin pertenecer a ninguna estructura virializada conteniendo galaxias 
donde la masa total pueda determinarse a traves de los metodos dinamicos habituales. Esa fue 
la posibilidad bajo la cual el modelo Einstein-de Sitter (J7 = 1) fue forzado a resguardarse de la 
evidencia observacional. Sin embargo, varios avances durante la ultima decada han confirmado un 
valor de O ~ 0.3. En los ciimulos de galaxias de mayor masa, la fraccion de la masa constituida 
por bariones deducida por la intensidad y el espectro de rayos X emitidos por el gas caliente es 
har - 0.06/1-3/2 (s^ ^ M_ White et al. 1993; D. A. White & Fabian 1994). Dada la densidad 
de bariones deducida de la teoria de nucleosmtesis primordial, ^^5 ~ 0.019/1^^ (p-e.. Buries & 
Tytler 1998), y el hecho de que la fraccion de masa barionica en los ciimulos mas masivos debe 
ser representativa de la fraccion media en el universo (White et al. 1993), deducimos un valor 
$7 ~ 0.3/i~^'2. Al mismo tiempo, las observaciones recientes de curvas de luz de supernovas tipo la 
deducen el mismo valor de fi, concluyendo ademas que la expansion del universo esta acelerando 
de la forma esperada en el modelo del universo con geometria espacial plana, donde la constante 
cosmologica proporciona la densidad de energia adicional necesaria para llegar a la densidad 
critica (Perlmutter et al. 1998, Riess et al. 1998). Otros metodos de medir el valor de O son, 
por lo general, consistentes con este resultado (por ejemplo, el valor $7 ~ 0.3 con una constante 
cosmologica A = 1 — i7 tambien es favorecido por el valor de la constante de Hubble y la edad de 
las estrellas mas antiguas en ciimulos globulares). Varias observaciones en el futuro proximo (tales 
como supernovas Tipo la y fluctuaciones en la radiacion cosmica de fondo) deberan permitir medir 
il con mayor precision, y clarificar si la reciente aceleracion del universo se debe a una constante 
cosmologica o al modelo mas general de un campo escalar con presion negativa (p.e., Peebles & 
Vilenkin 1998 y referencias incluidas). 

Una vez adoptamos el modelo cosmologico deducido a partir de esas observaciones, las 
predicciones de la teria MIF se ajustan bien a los datos observacionales sobre estructura a gran 
escala, tales como la funcion de correlacion espacial de las galaxias, la abundancia de ciimulos de 
galaxias, la evolucion en CR de esas cantidades, y las fluctuaciones en la radiacion cosmica de 
fondo recientemente detectadas. Para obtener esas predicciones, es preciso calcular el espectro de 
potencias de las fluctuaciones primordiales. 

El espectro de potencias de las fluctuaciones de densidad se obtiene en la teoria MIF 
calculando la evolucion en el regimen lineal de esas fluctuaciones una vez cruzan el horizonte 
(para articulos de revision, vease Efstathiou 1990, Bond 1996, Bertschinger 1996). En el limite de 
pequeiias escalas, las fluctuaciones entraron en el horizonte durante la epoca en que la densidad 



media del universe era dominada por radiacion. Los bariones estan entonces ligados a la radiacion, 
y la presion radiativa resulta en una oscilacion de las fluctuaciones, impidiendo el crecimiento 
gravitatorio. La materia invisible no interactiia con la radiacion, y por lo tanto aumenta sus 
fluctuaciones de densidad. No obstante, ese crecimiento gravitatorio es muy lento incluso en la 
ausencia de acoplamiento con la radiacion, cuando la radiacion domina la densidad de energia. 

Dado que la densidad de radiacion en el presente es 0,rad = (^T^ad/ (Pcritc^) = 2.47 x lO^^/i^^ 
(donde T^-ad = 2.76 K es la temperatura de la radiacion de fondo), la epoca de igualacion de las 
densidades de materia y radiacion corresponde al CR 1 + z = 4.04 x lO'^il/i^. La longitud del 
horizonte en esta epoca era del orden cH^^{z) = cHq il^^'^ (1 + z)^'^'^, y por lo tanto la longitud 
comovil era Lig = cHq 17^^"(1 + Zig)^^''^ ~ I5{i^h)^^ h^^ Mpc (donde utilizamos la notacion 
habitual para la constante de Hubble en el presente, Hq = 100/ikms~ Mpc~^, y H{z) es la 
constante de Hubble a CR z). A escalas comoviles mucho menores que el horizonte en la epoca de 
igualacion, la amplitud de fluctuaciones de densidad debe conservar la forma de las fluctuaciones 
primordiales cuando emergen del horizonte, ya que las fluctuaciones empezaron a crecer solamente 
a partir de la epoca de igualacion, y el crecimiento en la epoca anterior fue muy lento. Eso implica 
que la amplitud de fluctuaciones tiende a una constante para el caso mas habitual de invariancia 
de escala, con n = 1 (donde el espectro primordial es P{k) oc fc"). Pero en el limite de escalas 
mucho mayores, L ^ Lig, las fluctuaciones crecen en el regimen lineal proporcionalmente al factor 
de escala, a oc (1 + z)~^, desde el momento en que entran en el horizonte; dado que la longitud 
comovil del horizonte es proporcional a (1 + z)^^'^, el factor total de crecimiento a una epoca fija 
debe ser proporcional a L~^. 

Mostramos en la Figura 1 la dispersion en la fluctuacion relativa de la masa total 
contenida dentro de una esfera de radio comovil R, cuando la evolucion lineal de fluctuaciones 
primordiales se extrapola hasta el presente. Conviene recordar aqui que la distribucion de 
probabilidad de esta masa es Gausiana (dada la suposicion de que el campo de densidades 
es Gausiano), con la dispersion mostrada en la Figura 1. La linea solida es para el modelo 
il = 0.3, A = 0.7, n = 1, Hq = 65kms~ Mpc~ , y hemos normalizado las fluctuaciones a 
as =< o"^ >i/2 (i? = 8/i~^ Mpc) = 0.9. La linea de rayas es para el modelo fi = 1, A = 0, 
Hq = 50kms^ Mpc^ , y erg = 0.55. Las curvas han sido calculadas con las formulas presentadas 
por Hu & Sugiyama (1996). Ademas del comportamiento esperado en el limite de pequehas y 
grandes escalas, vemos tambien que el modelo de la constante cosmologica dispone de mayor 
potencia a grandes escalas, debido al mayor valor de la longitud caracteristica Lig oc (O/i)"-*^. 

La normalizacion del espectro de fluctuaciones, dada por el parametro cxg, debe considerarse 
como un parametro ajustable de la teoria. Su valor ha sido escogido aqui para reproducir 
aproximadamente la abundancia de cumulos de galaxias en el presente (Eke, Cole, & Frenk 1996). 
Desde que las fluctuaciones de intensidad de la radiacion de fondo de microondas fueron detectadas 
por el satelite COBE, la normalizacion viene fijada independientemente por esas fluctuaciones 
(Bennett et al. 1996; Bunn & White 1997). Esas dos medidas de la amplitud concuerdan para 
el modelo con constante cosmologica con il ~ 0.3, mientras que son claramente conflictivas si 



se supone ri = 1. Es notable el hecho de que cuando se toma el valor de $7 deducido de las 
observaciones de la geometria global del universo y la densidad promedio de materia, la teoria MIF 
predice correctamente la relacion entre las fluctuaciones de la radiacion de fondo y la abundancia 
de ciimulos de galaxias, y concuerda ademas con la correlacion espacial de galaxias (Jing, Mo, & 
Borner 1998 y referencias incluidas), asi como la evolucion de la abundancia de cumulos en CR 
(Bahcall & Fan 1998). 



3. Colapso gravitatorio no-lineal: formacion de halos 

En la teoria MIF, la materia invisible empieza a formar objetos virializados (o "halos") por 
colapso gravitatorio a partir de las escalas mas pequefias, cuando la fluctuacion de densidad 
alcanza un valor ~ 1, entrando en el regimen no lineal. Luego, cuando las fluctuaciones a escalas 
mayores llegan al colapso, halos de mayor masa se forman por mersion de halos de masa menor 
colapsados anteriormente. Solo podemos presentar aqui un resumen muy breve de la teoria de 
formacion de halos y galaxias; para una discusion mucho mas extensa, vease White (1996). 

Un modelo analitico de gran utilidad para entender el proceso del colapso no lineal, llevando 
a la formacion y mersion de halos, es el modelo de Press-Schechter (Press &: Schechter 1974). El 
modelo se basa en la solucion analitica del colapso de una perturbacion de densidad constante 
con simetria esferica, en cuyo caso el colapso al punto central sucede cuando se alcanza el valor 
de la sobredensidad extrapolada linealmente 6 = 6c = 3/5(37r/2)^/^ = 1.686 (Peebles 1980). 
La distribucion de la densidad promediada sobre una esfera de radio R es una Gausiana con 
dispersion (t{R) (Figura 1), y el modelo Press-Schechter consiste en suponer que la fraccion de 
la masa del universo que ha colapsado en halos de masa mayor que M = (47r/3)pi?^ es igual 
a erf{(5c/[2^'^cj(i?)]}. Notese que esa fraccion es en realidad dos veces la fraccion del volumen 
con sobredensidad superior a 5c en las condiciones iniciales. Aunque la introduccion del factor 
2 puede justificarse mejor con un tratamiento en el espacio de Fourier (Bond et al. 1991), el 
modelo de Press-Schechter debe considerarse linicamente como una aproximacion que se utiliza 
frecuentemente debido a su gran simplicidad, y especialmente al hecho de que sus predicciones se 
ajustan a resultados numericos de la abundancia de halos con bastante precision (p.e., Governato 
et al. 1998; a la practica, el parametro 5c puede considerarse ajustable a los resultados numericos, 
en vez de tomar el valor requerido para el colapso esferico). 

Una vez obtenida la abundancia de halos de masa mayor que M, se obtiene por diferenciacion 
la abundancia diferencial, y asignamos a cada escala R una dispersion de velocidades para el halo 
obtenida tambien por el modelo esferico: a^ = (37r/2)^'^ [H{zc)R\'^ 12, donde H{zc) es la constante 
de Hubble en el momento del colapso. Basicamente, esa relacion nos dice que la dispersion de 
velocidad de un ciimulo es del orden de la velocidad de expansion de Hubble sobre la region comovil 
desde la cual el ciimulo ha colapsado, en el momento del colapso. La dispersion de velocidad de 
un halo (o, analogamente, el cociente de la masa sobre el radio) es la propiedad de un halo que 
determina las cantidades observables en galaxias y cumulos, que son la dispersion de velocidad 



de las galaxias, la temperatura del gas difuso en equilibrio dinamico en el halo (determinada con 
espectros de rayos X), o la defleccion gravitatoria de luz de fondo (medida en lentes gravitatorias). 

En las Figuras 2(a,b), las tres lineas solidas gruesas nos dan la dispersion de velocidades de 
un halo formado a partir de una fluctuacion de densidad inicial igual a (1, 2, 3) — a, en funcion 
del CR, para los misnios dos modelos utilizados en la Figura 1 (para el modelo con constante 
cosmologica en la Fig. 2b, los valores de 5c y de la dispersion de velocidades correspondiente a 
cada escala R deben ser modificados debido al efecto de repulsion de la constante cosmologica, 
y el ritmo de crecimiento lineal de fluctuaciones es tambien distinto; vease p.e., Viana & Liddle 
1996). De acuerdo con la distribucion Gausiana, la fraccion de niasa que ha colapsado en halos de 
mayor dispersion de velocidades es, respectivamente, (64%, 10%, 0.6%). Tomemos, por ejemplo, el 
momento presente {z = 0) para el modelo $7 = 0.3 (Figura 2b). Ese modelo nos predice que un 
0.6% de la masa en el presente esta en objetos colapsados con dispersion de velocidad mayor que 
1000 km s^ . Esos halos se corresponden evidentemente con los ciimulos de galaxias mas masivos, 
y se formaron a partir del colapso de raras fluctuaciones a escalas grandes, i? ~ 15 h^^ Mpc, donde 
3a{R) = 5c- Las li'neas de rayas nos dan el valor de la masa de los halos (con un factor 10 en masa 
separando lineas sucesivas), e indican que la masa de uno de esos ciimulos es ~ 2 x IO^^Mq. La 
mayor parte de la masa del universo debe encontrarse en objetos formados a partir de fluctuaciones 
mas habituales: una fluctuacion \ — a produce un halo con dispersion de velocidades ~ 200kms~ , 
un valor tipico de los pequefios grupos de galaxias, tales como el grupo Local, donde residen la 
mayoria de las galaxias en el universo. La masa indicada en la Figura 2 se refiere siempre a la 
masa total del halo en el momento de su formacion. Por ejemplo, la galaxia de la Via Lactea pudo 
formarse a 2 ~ 2, en un halo con a ~ 150 km s^ de masa M ~ lO^^M©, y en el presente el Grupo 
Local esta colapsando en un halo mayor. Los dos modelos en la Figura 2 han sido normalizados 
para ajustar sus predicciones a la abundancia observada de ciimulos. 

La dispersion de velocidades puede relacionarse con la temperatura del gas difuso en el 
halo cuando esta en equilibrio hidorstatico: kT/ li = a^. Mostramos tambien la temperatura en 
el eje vertical. Hemos utilizado aqui la masa media para el caso de gas primordial ionizado: 
;U~0.6mH ~ 10~24g. 

Una propiedad general de la teoria MIF es que, a medida que las fluctuaciones colapsan a 
escalas cada vez mayores, los halos deben mergerse y formar nuevos halos con mayor dispersion 
de velocidad; ademas, el ritmo de crecimiento de la masa y la dispersion de velocidad es mucho 
mas rapido a alto CR que en el presente. Eso es una consecuencia directa de la forma del espectro 
de potencias en la Figura 1: a pequehas escalas, la amplitud de fluctuaciones es practicamente 
constante, y por lo tanto la masa de los halos aumenta rapidamente cuando esas fluctuaciones 
empiezan a colapsar. Pero a mayores escalas, la amplitud disminuye con mucha mas rapidez con 
la escala, y por lo tanto las fluctuaciones a escalas grandes colapsan mucho mas tarde. 



4. Formacion de galaxias 

Una pregunta central en la teoria de formacion de galaxias ha sido la siguiente: cual es la 
causa de la diferencia entre una galaxia y un ciimulo de galaxias? Evidenteniente, los cumulos de 
galaxias han colapsado gravitatoriamente, pero deducimos a partir de observaciones de rayos X 
que la mayor parte de la materia barionica esta en forma de gas difuso y caliente distribuido en 
el halo. La masa total del gas difuso en un ciimulo masivo puede superar las IO^^Mq, pero no se 
forman nunca galaxias de masa tan grande. 

Esa cuestion fue investigada por Binney (1977), Rees & Ostriker (1977), Silk (1977), y 
especialmente en el contexto de la teoria MIF por White & Rees (1978). En el colapso gravitatorio 
de un halo, la materia barionica se calienta inicialmente en ondas de choque, hasta la temperatura 
necesaria para permanecer en equilibrio hidrostatico. Luego, para poder condensar a una galaxia 
central, es precise que el gas pierda su energia termica por procesos de radiacion. Si el tiempo 
necesario para que toda la energia termica de los bariones sea radiada es corto comparado con 
el tiempo entre mersiones sucesivas de halos, entonces el gas va a poder concentrarse hacia el 
centre del halo, y el proceso de colapso gravitatorio puede proseguir hasta llegar a la formacion 
de estrellas (despues de fragmentacion de nubes de gas, y posiblemente de haber formado un 
disco debido a la conservacion del momento angular; vease Fall & Efstathiou 1980, Fall & Rees 
1985) o de un objeto masivo central, dando lugar a un niicleo activo. En cambio, si el tiempo 
de enfriamiento es demasiado largo, nuevas mersiones para formar halos de mayor masa van a 
calentar de nuevo el gas, manteniendolo en estado difuso. 

La condicion que debemos analizar es, por consiguiente, que el tiempo de enfriamiento del gas 
sea igual a la edad del universo (que es aproximadamente el tiempo entre mersiones sucesivas). 
Los procesos de enfriamiento mas importantes son la radiacion bremsstrahlung, y la excitacion y 
ionizacion colisional de iones. A temperaturas T ^ lO^K, el enfriamiento por hidrogeno y helio 
domina, pero en el rango lO^K < T < 10^ K los iones de elementos pesados son importantes 
(p.e., Gaetz & Salpeter 1983). En las figuras 2(a,b), las lineas de puntos indican la dispersion de 
velocidad de halos en los que el tiempo de enfriamiento es igual a la edad del universo, para tres 
valores de la metalicidad del gas, (0,0.1, 1)Zq. Por encima de las lineas de puntos, el tiempo de 
enfriamiento es demasiado largo y se predice que el gas caliente se mantiene en el halo, mientras 
que por debajo de las curvas el gas puede enfriarse y formar una galaxia. El ritmo de enfriamiento 
por radiacion se ha calculado suponiendo que una fraccion i^h/Q de la masa del halo esta en forma 
de gas difuso, donde Qh = 0.019/i~^, de acuerdo con las mediciones de la abundancia de deuterio 
(Buries &: Tytler 1998), y suponiendo el valor de la sobredensidad media de los halos en el modelo 
del colapso esferico (igual a 187r^ para fi = 1). Si las galaxias pueden formarse en halos por debajo 
de las lineas de puntos, obtenemos la prediccion de que hasta z ;^ 4, el colapso de un halo debio 
resultar siempre en la formacion de una galaxia central, pero a menores CR los halos mas masivos 
pueden contener solamente galaxias formadas anteriormente, y la mayor parte del gas difuso no 
puede enfriarse. Mencionamos aqui que, incluso en los cumulos de mayor masa, una parte del gas 
en la region central (con mayor densidad) puede radiar mas rapidamente; la razon por la cual eso 



no da resultado a grandes tasas de formacion estelar en las galaxias centrales de cuniulos masivos 
es una cuestion que pernianece sin resolverse (Fabian, Nulsen, & Canizares 1994). 

Vemos, pues, que la fisica de enfriamiento del gas proporciona una explicacion satisfactoria 
del limite superior de la masa de las galaxias. Otra cuestion distinta es el limite inferior, y eso nos 
lleva a una transicion fundamental en el universo para la historia de formacion de las galaxias: la 
reionizacion del medio intergalactico. 



5. La reionizacion del universo 

En la teoria de la Gran Explosion, la materia del universo debe formar atomos por vez 
primera cuando la temperatura de la radiacion de fondo disminuye hasta T ~ 3000 K, a z ~ 10^. 
La materia intergalactica permanece posteriormente en estado atomico, hasta el momento en que 
el colapso no-lineal de las perturbaciones lleva a la formacion de objetos que emiten radiacion 
ionizante, y pueden ionizar de nuevo el gas difuso. 

Las observaciones del espectro de fuentes luminosas a alto CR (generalmente, cuasares) 
demuestran que el medio intergalactico fue ionizado anteriormente a z ~ 5. La luz de una fuente 
a longitudes de onda menores que la linea Lya de hidrogeno atomico puede ser absorbida por 
un atomo a lo largo de su trayectoria a traves del universo, en el punto en que la longitud de 
onda se ha corrido hasta coincidir con la linea Lya . Si una parte importante de la densidad 
barionica media del universo estuviera distribuida por el espacio en forma atomica, la profundidad 
optica de absorcion seria enorme {^ 10^ a z = 3), con lo cual el flujo deberia disminuir a cero 
abruptamente a la longitud de onda de Lya (Gunn & Peterson 1965). Lo que se observa en 
realidad es la presencia de multiples lineas de absorcion que producen un decremento neto del flujo 
de solamente un ~ 30% a z = 3. Eso implica que el medio intergalactico esta altamente ionizado, 
con una fraccion neutra de ~ 10^^. Esa diminuta fraccion del hidrogeno intergalactico puede 
explicar satisfactoriamente las propiedades de las lineas de absorcion observadas (denominadas 
usualmente como el "bosque Lya "), originadas en las variaciones de la densidad del gas que no 
ha colapsado todavia en halos virializados de gran densidad. Varios trabajos recientes, utilizando 
modelos analiticos y simulaciones numericas, han mostrado que la teoria MIF predice de forma 
generica que el medio intergalactico ionizado debe dar lugar a ese bosque de lineas (vease McGill 
1990; Bi, Borner, & Chu 1992; Bi 1993; Miralda-Escude & Rees 1993; Gen et al. 1994; Zhang et 
al. 1995, 1998; Hernquist et al. 1996; Miralda-Escude et al. 1996). El hecho de que la densidad 
del gas neutro es proporcional al cuadrado de la densidad del gas cuando se establece el equilibrio 
entre recombinacion y fotoionizacion por el fondo cosmico ionizante es la causa de que un medio 
continuo da lugar a un espectro de absorcion con la apariencia de "lineas" individuales, cada vez 
que se encuentra un maximo de la densidad del gas a lo largo de la linea de visual. El espectro 
de absorcion tiene una escala de suavizamiento natural de ~ 20kms^ , debido a la dispersion de 
velocidad termica del gas fotoionizado (con temperatura T ~ 2 x 10'^ K). El tamaiio transversal de 
las estructuras de absorcion Lya debe ser del orden de su dispersion de velocidad multiplicada por 



la edad del universo, ya que tales estructuras no han tenido todavia tiempo de colapsar y de llegar 
a un equilibrio hidrostatico; observaciones recientes del tamano transversal en pares de cuasares 
(Bechtold et al. 1994; Dinshaw et al. 1994) confirman esa prediccion. 

Que causo la reionizacion del universo? Dos mecanismos pueden ionizar el gas intergalactico: 
fotoionizacion, o ionizacion colisional una vez el gas ha sido calentado por ondas de choque 
provenientes de alguna explosion de gran energia. La fotoionizacion es el metodo mas eficiente 
para ionizar el medio de baja densidad, ya que requiere menos energia, y la radiacion se transporta 
con gran eficacia a todas las regiones del espacio. Los dos tipos de fuentes de radiacion ionizante 
conocidas que pueden formarse tan pronto como los primeros halos colapsan a pequeiias escalas 
son estrellas y micleos activos (producidos por agujeros negros masivos en el centre de una 
galaxia). La cantidad de estrellas necesaria para ionizar el universo entero es solamente una parte 
muy pequeria de todas las estrellas que se han formado hasta el presente, y puede relacionarse 
facilmente con la metalicidad media producida. Ya que tanto los elementos pesados provenientes 
de explosiones supernova como la radiacion ionizante son producidas por estrellas masivas, la 
razon de esas dos cantidades esta relativamente fijada: caracteristicamente, una estrella de 30Mq 
fusiona unas 4Mq de hidrogeno en la secuencia principal, obteniendo una energia ~ O.OSMqc^, 
de la cual ~ O.OIMqc^ se emite en fotones ionizantes. Esa misma estrella produce unas 5Mq 
de elementos pesados en la explosion supernova. Si los fotones son absorbidos por hidrogeno 
intergalactico (necesitando una energia de ~ 20 eV para cada ionizacion, o una fraccion 2 x 10~* 
de la masa-energia en reposo), la masa ionizada es de ~ 5 x 10^ Mq, y por lo tanto la metalicidad 
media aumenta solo hasta 10~^ una vez se ha emitido un foton ionizante para cada barion. Dado 
que la metalicidad media en el universo presente es mucho mayor, es evidente que la primera 
generacion de estrellas puede facilmente ionizar todo el universo (p.e., Couchman & Rees 1986). 

Evidentemente, los cuasares pueden ser tambien las fuentes dominantes para la reionizacion, 
dado que la eficiencia en convertir la masa acretada por un agujero negro en radiacion ionizante 
es generalmente mucho mayor que para estrellas. Las observaciones de la abundancia de cuasares 
y de la intensidad del fondo cosmico de radiacion ionizante a z ^ 3 indican que los cuasares son 
probablemente los mayores contribuyentes de esa radiacion cosmica, pero la naturaleza de las 
fuentes causantes de la reionizacion a mayor CR es todavia incierta (Miralda-Escude & Ostriker 
1990; Madau 1991, 1992, 1999; Haardt & Madau 1996; Ranch et al. 1997; Haiman & Loeb 1997, 
1998; Miralda-Escude, Haehnelt, & Rees 1999). 



6. Las primeras galaxias 

Las primeras galaxias donde pudieron formarse las primeras estrellas y cuasares surgieron, 
en la teoria MIF, en los primeros halos que colapsaron donde el gas pudo enfriarse. El hidrogeno 
atomico solamente empieza a radiar en lineas de excitacion cuando la temperatura supera los 10 
K; a temperaturas inferiores, las colisiones con electrones termicos no tienen nunca la energia 
necesaria para excitacion a los niveles atomicos n = 2. El hidrogeno molecular proporciona la 
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linica fuente de enfriamiento a nienores temperaturas. En la materia primordial, en ausencia total 
de elementos pesados, el hidrogeno molecular se forma a partir del ion H^ (Saslaw & Zipoy 1967), 
el cual resulta a su vez de colisiones de hidrogeno con los electrones del residuo de ionizacion 
que permanece despues de la epoca de recombinacion (Peebles 1968). Ese proceso de formacion 
de moleculas es muy ineficaz (en el presente, el hidrogeno molecular se forma en la superficie 
de granos de polvo interestelar) , y solo una pequeha fraccion del hidrogeno forma moleculas. 
El enfriamiento resultante solo es suficiente para resultar en la formacion de galaxias cuando la 
temperatura de un halo supera los ^ 2000-ftr (Tegmark et al. 1997; Abel et al. 1998). 

En las Figuras 2(a,b), la li'nea solida a mayor CR indica la minima temperatura para que 
un halo formado a CR z pueda enfriar su gas y formar una galaxia, que hemos reproducido de 
Tegmark et al. (1997). Vemos que la teoria MIF, con los parametros obtenidos de observaciones 
mencionadas antes, predice que las primeras estrellas pudieron formarse a z ~ 20, en halos de 
dispersion de velocidades de -^ 5kms~ y masa total ~ 10 Mq. Esa primera generacion de 
galaxias, formadas a traves del enfriamiento por hidrogeno molecular, formo probablemente una 
cantidad muy pequeha de estrellas debido a la ineficacia de esta forma de enfriamiento (vease 
tambien Haiman, Rees, & Loeb 1997). En esta epoca, el ritmo de mersiones es muy elevado y 
la masa y dispersion de velocidad de los halos aumenta rapidamente con el tiempo; asi pues, 
inmediatamente despues de la formacion de esas primeras estrellas, el enfriamiento atomico en 
halos de ~ W^Mq empieza a ser activo, dando lugar a una nueva generacion de galaxias que debio 
reionizar el universo. 

Una vez el medio intergalactico es reionizado, el enfriamiento es suprimido por la 
fotoionizacion, que calienta el gas y disminuye la abundancia de atomos que pueden ser excitados 
colisionalmente. La linea solida delgada a bajo CR nos da la temperatura en la cual el tiempo 
de enfriamiento es menor que el tiempo de Hubble para el gas ionizado (hemos adoptado aqui 
el modelo de Haardt & Madau 1996 para el fondo de radiacion ionizante). Las galaxias pueden 
formarse solamente por encima de esta linea una vez la reionizacion se ha completado. En realidad, 
la fotoionizacion tiene tambien otro efecto: el calentamiento del gas durante la reionizacion puede 
producirse antes del colapso, y luego la temperatura del gas sube adiabaticamente cuando la 
densidad aumenta si el enfriamiento no es importante (recordemos que las lineas indicando un 
tiempo de enfriamiento igual a la edad del universo son para una sobredensidad de 187r^ ~ 178, 
el valor obtenido en el modelo esferico en el momento de la virializacion, y que el tiempo de 
enfriamiento es mayor a menor densidad). Debido a este efecto, el gas en halos con a ^ 30kms~ 
no puede en general enfriarse para formar galaxias (p.e., Thoul & Weinberg 1996). Otro efecto 
que puede disminuir la eficiencia de formacion de galaxias en halos de baja dispersion de velocidad 
es el hecho de que la energia liberada por el propio proceso de formacion estelar puede crear un 
viento galactico con energia suficiente para calentar y expulsar el gas en el proceso de acrecion 
(p.e., Dekel & Silk 1986). Ese proceso parece ser necesario para evitar un exceso de galaxias de 
baja luminosidad comparado con las observaciones (White & Frenk 1991; Navarro & Steinmetz 
1997). 
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7. Observaciones de galaxias a alto CR 

Grandes avances en cosmologia observacional han tenido lugar recientemente, con el 
descubrimiento de gran niimero de galaxias a alto CR. Una de las tecnicas de mayor exito es la 
seleccion fotometrica de objetos debiles para deteccion de la brecha de Lyman ("Lyman break"; 
Guhathakurta, Tyson, &: Majewski 1990; Steidel et al. 1996). Todas las galaxias, cuya luz resulta 
de la superposicion de los espectros de muchas estrellas, muestran en su espectro una cai'da abrupta 
del flujo a la longitud de onda del limite de Lyman, como vemos en los modelos de espectros de 
distintas edades con ritmo constante de formacion estelar en la Figura 3, que reproducimos de 
Bruzual & Chariot (1993). La brecha se produce en la atmosfera de las estrellas masivas, que 
producen la mayor parte de la radiacion ultravioleta. Ademas, la presencia de hidrogeno atomico 
en el medio interestelar de la galaxia resulta generalmente en la absorcion de la mayor parte de 
los fotones ionizantes, aumentando la amplitud de la brecha. A alto CR, la brecha de Lyman se 
corre hasta longitudes de onda en el visible, lo que permite seleccionar las galaxias a partir de sus 
colores. Por ejemplo, un objeto azul en V-R y muy rojo en B-V debe tener la brecha entre las 
bandas B y V, situandolo a z ~ 4. 

El mimero de galaxias detectadas por este metodo, de magnitudes en el rango 23 ^ B ^ 29, 
ha permitido la primera medicion de la funcion de luminosidad y de la tasa global de formacion 
estelar a alto CR, que resulta ser mucho mayor que la actual (Madau et al. 1996; Steidel et al. 
1999). La cantidad total de estrellas formadas en estas galaxias, en el intervalo 1 ^ ^ ^ 5, puede 
dar cuenta de la mayor parte de la poblacion estelar vieja en el universo presente. Sin embargo, 
esas estimaciones estan sujetas todavia a varias incertidumbres: en general, solamente la tasa de 
formacion de estrellas masivas se puede deducir de esas observaciones, y existe la posibilidad de 
que la funcion de masa inicial fuera distinta a alto CR. Por otra parte, la tasa total de formacion 
estelar podria ser mucho mayor si galaxias de baja luminosidad que permanecen por debajo de los 
limites de deteccion dominaran la emision total, o si la mayor parte de la radiacion ultravioleta es 
absorbida por polvo interestelar y reemitida en el infrarojo lejano (p.e., Calzetti 1999 y referencias 
incluidas) . 

Efectivamente, la radiacion reemitida por polvo constituye otra forma de deteccion de galaxias 
a alto CR. Recientemente, la radiacion de fondo en el infrarojo lejano, debida a la combinacion de 
todas las galaxias, ha sido detectada por COBE (Fixsen et al. 1998; Hauser et al. 1998), y fuentes 
individuales se han descubierto con el nuevo detector SCUBA (Eales et al. 1998; Hughes et al. 
1998; Small et al. 1997), sugeriendo que gran parte de la radiacion estelar a alto CR pudo ser 
procesada por polvo. 

El origen de esas galaxias a alto CR puede comprenderse simplemente en la teoria MIF a 
partir de la Figura 4. Las li'neas punteadas nos indican aqui la magnitud aparente de una galaxia 
(en el sistema AB, donde la magnitud es AB = —48.6 — 2.5logio{fu), y el flujo fi, se expresa en 
unidades cgs), formada en un halo de dispersion de velocidad a en funcion de z, para un modelo 
de maxima luminosidad en que todos los bariones contenidos en el halo forman estrellas en un 
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intervalo de tiempo igual a la mitad de la edad del universe en el moniento del colapso. En otras 
palabras, el modelo representa la maxima eficiencia de formacion estelar posible, donde toda la 
materia forma estrellas en un tiempo del orden del tiempo dinamico del halo de materia invisible 
(vease Miralda-Escude & Rees 1998 para mas detalles del modelo). Las galaxias de mayor masa 
pueden empezar a formarse a z ~ 4, cuando el gas puede enfriarse segiin la figura 2b. El flujo 
maximo de esas galaxias corresponde a AB ~ 22. La mayor parte de las galaxias de la brecha de 
Lyman son algo mas debiles, como es de esperar cuando el ritmo de conversion de gas a estrellas 
es menos eficiente, y cuando una parte de la radiacion ultravioleta es absorbida por polvo. Vemos 
tambien que las galaxias mas masivas pueden formarse mas facilmente cuando el gas ha sido ya 
enriquecido, debido al mayor ritmo de enfriamiento a alta metalicidad. 

Puesto que las galaxias detectadas a alto CR de mayor luminosidad estan asociadas con las 
mayores fluctuaciones de densidad a las mayores escalas de colapso gravitatorio, su correlacion 
espacial deberia ser mucho mayor que la correlacion de la masa, como se espera de los altos picos 
de densidad en un campo Gausiano (Kaiser 1984). Esa correlacion ha sido detectada, y es en 
general consistente con las expectativas en la teoria MIF (Adelberger et al. 1998; Giavalisco et al. 
1998; Kauffmann, Nusser, & Steinmetz 1997). Los estudios de la correlacion de galaxias a alto CR 
abren un nuevo abanico enorme de posibilidades para observar la evolucion de estructura a gran 
escala, que hemos solamente empezado a investigar. 

A mayor z, la absorcion del flujo ultravioleta a A < 1216 A debido al bosque Lya aumenta 
rapidamente, de tal forma que a z ^ 5, la caida de flujo en Lya se convierte en la caracteristica mas 
importante para seleccionar objetos a este CR. Asi pues, la tecnica para encontrar galaxias a z > 5 
va a ser muy parecida a la de la brecha de Lyman, sustituyendo esta por el hoyo de Gunn-Peterson 
(Gunn & Peterson 1965). Aunque inicialmente el termino del hoyo de Gunn-Peterson se utilize 
unicamente para referirse a la absorcion en Lya producida por el medio intergalactico atomico 
antes de la reionizacion, el medio ionizado puede dar lugar tambien a un hoyo cuando la fraccion 
neutra es suficientemente grande para que todas las partes del medio intergalactico, incluso las 
de menor densidad en los vacios, absorban esencialmente todo el flujo. Debido al aumento de la 
densidad media del gas en el universo, y al mayor ritmo de recombinacion a alto CR, es de esperar 
que practicamente todo el flujo a longitudes de onda menores que la linea Lya sea absorbido a 
z '^6, incluso si la reionizacion ocurrio a mayor CR (Miralda-Escude et al. 1999). 

Finalmente, otra tecnica importante de deteccion de galaxias a alto CR es mediante la linea 
de emision de Lya (Thompson et al. 1995; Thommes et al. 1998; Meisenheimer et al. 1998; Hu, 
Cowie, (fe McMahon 1998). En regiones de formacion estelar, la mayor parte de la radiacion 
ionizante emitida por estrellas jovenes es generalmente absorbida por hidrogeno interestelar, y la 
energia se reemite en fotones de recombinacion, siendo Lya la linea mas brillante. La biisqueda 
de lineas de emision de objetos a alto CR tiene la ventaja de que el fondo de cielo puede reducirse 
observando solamente en una banda estrecha de longitudes de onda, especialmente cuando se 
seleccionan longitudes de onda de buena transparencia atmosferica (especialmente importante a 
2^5, cuando la linea Lya se corre hasta el infrarojo). 
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Es evidente a partir de las Figuras 2 y 4 que la dificultad en la deteccion de galaxias a CR 
progresivamente mayor aumenta rapidamente para ^; ^ 5, puesto que el flujo de las galaxias debe 
disminuir no solo debido al mayor CR, sino a la baja luminosidad de las primeras galaxias, debida 
a su menor masa comparado con las galaxias actuales. Esa prediccion de la teoria MIF sugiere 
que en el future, las fuentes detectadas a mayor CR podrian ser supernovas, las cuales debieron 
ocurrir tan pronto como las primeras estrellas se formaron (Miralda-Escude & Rees 1997). La 
proyectada mision NGST podri'a detectar supernovas hasta z ~ 10 (Stockman &: Mather 1997). 
Una posibilidad interesante que podria permitir acelerar el estudio observacional de las primeras 
estrellas es que los estallidos de rayos gamma, con sus brillantes contrapartidas opticas (p.e., 
Metzger et al. 1997), sean un fenomeno asociado con estrellas masivas, ocurriendo por consiguiente 
a los mayores CR donde existian estrellas. Las contrapartidas opticas serian detectables en 
observatories terrestres, corridas al infrarojo, mientras que en el visual no habria contrapartida 
debido al hoyo de Gunn-Peterson y absorcion por fotoionizacion. En el caso de los estallidos de 
rayos gamma, el reto observacional podria consistir en identificar una pequeiia fraccion de esos 
eventos a mayor CR que cualquier otra fuente conocida, y distinguirlos de estallidos a menor CR 
pero con un gran enrojecimiento debido a polvo interestelar cerca de la fuente, que pueden tener 
caracteri'sticas fotometricas similar es. 



8. Conclusiones 

La cosmologia observacional esta entrando en una nueva etapa de descubrimiento, con nuevas 
tecnicas para la deteccion de galaxias debiles que empujan la frontera de alto CR hacia la epoca de 
la formacion de las primeras galaxias. Al mismo tiempo, la medicion precisa de las fluctuaciones en 
la radiacion de fondo, la realizacion de nuevos escrutamientos de CR de galaxias, y la continuacion 
de la busqueda de supernovas a alto CR para la medicion de la geometria cosmica, prometen fijar 
los parametros del universo y del modelo cosmologico. Hemos visto que la teoria MIF reproduce 
con gran exito las observaciones acumuladas hasta el presente sobre estructura a gran escala. 
Sin embargo, el estado actual de la teoria deja muchas preguntas abiertas: cual es la naturaleza 
de la materia invisible? Existe realmente una "energia de vacio" que da cuenta de la densidad 
necesaria para alcanzar la densidad critica? Cual es la naturaleza de esta energia de vacio, cual es 
su ecuacion de estado, y por que existe? Que proceso genero las fluctuaciones de densidad? Que 
determino su amplitud? Son las fluctuaciones primordiales exactamente adiabaticas y Gausianas, 
y es su espectro perfectamente invariante de escala, o existen pequeiias desviaciones de esta simple 
hipotesis? Esas preguntas, que nos llevan al misterio de la epoca de inflacion y al origen del 
universo, ocuparan probablemente el centro de interes en el futuro de la cosmologia observacional, 
cuyo avance permitira tambien descifrar la historia de la formacion de galaxias, desde el colapso 
de las primeras estrellas hasta el presente. 

Quisiera agradecer David Weinberg por proporcionar un codigo para calcular el espectro de 
potencias, y Gustavo Bruzual, Stephane Chariot y Max Tegmark por permitir la reproduccion de 
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Haehnelt, Martin Rees y David Weinberg. 
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Fig. 1. — Dispersion en la sobredensidad media dentro de una esfera de radio R en el campo lineal 
de densidad para la teoria MIF con constante cosmologica, con dos valores de Q. 
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Fig. 2a. — Las tres Imeas solidas gruesas dan la dispersion de velocidades de un halo colapsado a 
CR z, provenientes de fluctuaciones {l,2,3)a en el modelo Press-Schechter. La masa total de los 
halos es indicada por la Imea de rayas. Las lineas de puntos indican la maxima masa para que el 
gas pueda enfriarse, para tres metalicidades distintas, y las lineas solidas finas indican la minima 
temperatura para enfriamiento atomico y molecular. El modelo cosmologico es 17 = 1, A = 0, 
h = 0.5. 
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Fig. 2b. — Igual que en la Figura 2a, para el modelo il = 0.3, A = 0.7, h = 0.65. 
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Fig. 3. — Modelo del espectro de una galaxia con ritmo constante de formacion estelar, con las 
edades en aiios indicadas. 
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Fig. 4. — Figura analoga a la figura 2, donde las lineas de puntos dan la magnitud AB en una 
banda ultravioleta suponiendo que toda la niasa barionica en el halo forma estrellas durante un 
tiempo igual a la mitad de la edad del universo cuando el halo colapsa. 



